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Abstract 

A strong local form of the "4/3-law" in turbulent flow has been proved recently by 
Duchon and Robert for a triple moment of velocity increments averaged over both a bounded 
spacetime region and separation vector directions, and for energy dissipation averaged over 
the same spacetime region. Under precisely stated hypotheses, the two are proved to be 
proportional, by a constant 4/3, and to appear as a nonnegative defect measure in the local 
energy balance of singular (distributional) solutions of the incompressible Euler equations. 
Here we prove that the energy defect measure can be represented also by a triple moment of 
purely longitudinal velocity increments and by a mixed moment with one longitudinal and 
two tranverse velocity increments. Thus, we prove that the traditional 4/5- and 4/15-laws of 
Kolmogorov hold in the same local sense as demonstrated for the 4/3-law by Duchon- Robert. 
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1 Introduction 

Recently, Duchon and Robert have estabUshed an energy balance relation for distributional 
solutions of the three-dimensional (3D) incompressible Euler equations. Their balance relation 
contains a "defect" or "anomaly" term, with an interesting connection to turbulence theory. 
Since the work of Duchon-Robert provides the point of departure of the present paper, it is 
appropriate to describe their theorems briefly here. For our purposes, there are three main 
results: 

First, if u € -^^^([0, T] x T'^) is a weak solution of the incompressible Euler equations on the 
3-torus T^, then it is proved in ^ (Proposition 2) that the following local balance holds in the 
sense of distributions: 

9*(^|up) + V.[(^|up +p)u] = -D{u). (1.1) 

Here D(u) is a defect distribution which for classical solutions vanishes identically, implying 
local energy conservation. Duchon and Robert also establish various expressions for the defect 
term. In particular, they have shown that 

D{u) = v-iim\ [ v^sii)-Su{e)\6u{e)\'^ d^e (1.2) 

where V— lim means limit in the sense of distributions on [0,r] x T^, with ip£C^(T^), even, 
nonnegative with unit integral, ^p^{x) = e~^ip{x/e), and 6\i{x,t;£) = u{x + £,t) — u(x,t). 
This expression is remarkable because it is closely connected with an exact result in turbulence 
theory, the so-called "Karman-Howarth-Monin relation": 

V-((5u(^)|<5u(£)|2) = -4e. (1.3) 

See section 6.2.1. In this relation, e = i^dVu^^l^) is the mean energy dissipation for a 
Navier-Stokes solution u^, which is assumed to remain finite as viscosity u tends to zero. We 
see that this result is essentially equivalent to the statement that the Duchon-Robert defect 
satisfy (D) = e. 
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For a general distributional solution of Euler, there need be no connection of the defect 
with viscous dissipation, nor need it even be true that D{u) > 0. However, a second theorem 
(Proposition 4) in Q states that, if u'^ is a Leray solution of incompressible Navier-Stokes 
equation for viscosity u and if u*^ — > u strong in L^([0, T] x T'^) as ^ 0, then 

D{u) = V- lim[i/| Vu'^p + D{m'')] (1.4) 

Since it is well-known that D(u'^) > for Leray solutions, thus also -D(u) is a nonnegative 
distribution, i.e. a Radon measure (Q, Example 12.5). This theorem may be paraphrased as 
saying that strong viscosity solutions of the incompressible Euler equations are also dissipative 
solutions. 

Finally, Duchon and Robert, under an additional hypothesis, establish an even simpler form 
of the defect distribution. With uj unit Haar measure on 5^, they define the function 

S{u,e){x,t) :=-[ duj{i)SuL{x,t;e)\Su{x,t;i)\'^. (1.5) 

in L^([0,T] X T^). Here 5uL{^,t;£) = £-6u{x,t;£) is the longitudinal velocity increment. As- 
suming that the following limit exists 

5(u)(x,f) := P-lim5(u,£)(x,t) (1.6) 

Duchon and Robert Q, Section 5 show that 

5(u) = -^D{u). (1.7) 

This is a rigorous form of another well-known relation in turbulence theory, sometimes called 
the "4/3-law": 

{6uLimn{£)\^) ~ -^ei. (1.8) 

See Q . It is well-known that the Karman-Howarth-Monin relation reduces to the 4 / 3-law under 
conditions of local isotropy. This is achieved here by the angle average over the sphere in ( |1.5[ ). 
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Such relations as ( [L.3| ) and ( |1.8| ) in turbulence theory go back to the original work of A. 
N. Kolmogorov H]. However, Kolmogorov in fact proved a relation involving only longitudinal 
velocity increments, the so-called "4/5-law": 

{[6uL{e)f) ~ -^e-i (1.9) 

This was established from the Navier-Stokes equations, under conditions of statistical homo- 
geneity and local isotropy and with the assumption that energy dissipation remains finite in the 
zero-viscosity limit. It is our purpose here to establish an expression for the Duchon-Robert 
energy dissipation anomaly -D(u) with exactly the form of Kolmogorov's law. Our proof yields 
another expression which is related to the so-called "4/15-law": 

{6ul{£)\Sut{£)\^) ^ -^ei (1.10) 
15 

Here Suxii) = t-6u{i) is a transverse velocity increment, with t any unit vector orthogonal 
to £. This relation is known to hold under the same conditions as the 4/5-law. The rigorous 
derivation of the 4/5- and 4/15-laws given here, under precisely stated assumptions, yields a 
result with a wider domain of validity than that of some previous rigorous derivations, such as 
that of Nie and Tanveer In particular, the form of the 4/5-law established here — like the 
Duchon-Robert version of the 4/3-law — is local, in the sense that it relates third-order moments 
of velocity increments and viscous dissipation averaged over the same bounded spacetime region, 
not necessarily large. 

In the following section we prove our main theorem. In a final discussion section we discuss 
briefly its physical significance and a compare it with previous results. Let us now make just 
a few remarks on notations: The symbol will be used below for L^([0,T] x T^). If F(x, t) 
is any spacetime distribution, we denote F^(x,t) = (99^ * F){x,t), where G C^{T^) and * 
is space convolution. Thus, remains, for each fixed x, a distribution in t. Often below, as 
above, we omit the variables (x, t) where their presence is clear from the context. 
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2 The Main Theorem 

We prove the following: 

Theorem 1 Let u G L^{[0,T] x T^) be a weak solution of the incompressible Euler equations 
on the 3-torus . Let (p{e) be any C°° function with compact support, nonnegative with unit 
integral, spherically symmetric, and let ip^{£) = e~^<f{£/e). Finally, define longitudinal and 
transverse velocity increments as 

6uLix,t;e) = (i«)£)5u(x,t;€), ()uT(x,t;£) = (1 - ^0 i)5u(x, t; £), (2.1) 

Then, the following functions in L^{[0,T] x T^) 

Diiu) = ^J^^ dH I V^^(£)-5u(£)|5ul(£)|2 + j^^{i)5uL{£)\6uT{£)\''^ , (2.2) 

and 

D't{u) = ^J^^ dH ^V^'{iy5n{£)\5uTm^ - jip' {t)6uL{muT{tt^ , (2.3) 

both converge in the sense of distributions as e ^ to D(u), where the latter is the defect 
distribution in the local energy balance for u, 

dt{\\n\') + V.[(i|u|2 +p)n] = -D{n), (2.4) 

established earlier by Duchon-Robert. 

We shall prove this theorem in several steps. 

The idea of the proof is to consider separate balance equations for the longitudinal and 
transverse components of the energy. We define first longitudinal and transverse velocities 
relative to a vector £: 

UL{x,t;£) = {i(^£)-u{^ + e,t), ur(x,i;^) = (1 -£(8)^)-u(x + £,t). (2.5) 

Of course, ul + = u. We define also mollified versions 

u5f(x,t)=/" d^i (p'{e)uxi^,t;£), X = L,T (2.6) 
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It is easy to see that these satisfy the equations 

9tu5^ + V-(u®uxr = -V-n5^, X = L,T (2.7) 
distributionally in time, where 

ni(x,t)=/ d^i ip''{e){i(s>e)p{yL+£,t), n|.(x,t)=/ d^i ip^{£){i-e(s>e)p{yL+£,t). {2.8) 

These equations can be simpHfied by the observation that 

V-n5,(x,t) = Vp5,(x,t), X = L,T (2.9) 
where p^-.Pj' are scalar functions defined by 

p'x{x,t)= [ ip'xii)pix + i,t), X = L,T (2.10) 

with 

iPL{£) = m-M(), M^) = ^l^°^^de'. (2.11) 

Note that ipi, ipx are compactly supported and (7°° everywhere except at 0, where they have a 
mild (logarithmic) singularity. To prove (|2.9| ), we use the elementary relation = 1 — £ (E> £■ 
A simple computation then gives, for example, 

V-Ill{^,t) = - j^^dH |^(£) + ^^-(£)|£p(x + £,t), (2.12) 

From its definition, 

Vv^L(^) = {^(^) + ^9^^(^)}i (2.13) 

This gives easily V-III^ = Vp|^ in the sense of distributions. Because + = p^l and 
p\ + Pj^ = p^, this yields also the relation V-II^^ = Vpj.. Finally, we obtain the simpler 
equations 

9tU^ + V-(u®ux)" = -Vp^, X = L,T (2.14) 

for and u|,. 
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We next observe that both and u|i are divergence-free. In fact, a computation hke that 
above shows that 

V-ui(x,i) = -^^^3^ |^(£) + ^^-(£)|£.u(x + £,t), 

= - [ d^£V(pl{£)-u{x + £,t). (2.15) 

If we integrate against any smooth test function ip{:>i.), we then get 

/ d^x V(x)V-ul(x,t) = / d^£ ^l{i)V-{'ip*u){i,t). (2.16) 

The latter is zero, since u is divergence- free (in the distributional sense). Thus, V-u^^ = 0. 
Since -|- Uj, = u^, then also V-Uj. = 0. 

From the equations ( p. 141 ) for X = L,T, the incompressibility conditions for and Uj,, 
and the Euler equations for u (in distribution sense), we derive the following balance equations 

2dt{u-ul) + V-[2(u.ui)u + {{ul-ul)uY - {ul-ulYu + 2pul + 2plu] = -^DUu), (2.17) 

and 



2dt{u-u^T) + V-[2(u-u|.)u + ((ut-ut)u)^ - (ut-ut)^u + 2puf. + 2p|.u] = — D^{u). (2.18) 

3 

The basic identities used to derive these equations are 

d^f lv^'{ey6u{£)\5nL{£)\^ + j^'{£)6ul{£)\Sut{£)\^ 

■ -^{hijV'{l)}5ui{£)5uj{£)5uk{£) 
d d 

-^[((uL-UL)ufc)^ - (uL-UL)^lifc] + 2Ui—^[{{uLiUkf - U%Uk], (2.19) 
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and 
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{ W^'{l)-5vi{£)\5viT{£)? - jH^'{t)5uL{£)\5viT{£)? 
■ ^''^ - hij)^'mSM£)6uj{£)6uk{£) 

J3 Oik 

d d 

-[((UT-Ur)ufc)^ - {uT-UTYuk] + 2Ui-^[{{uTiUky - UTiUk]. (2.20) 



dxk dxk 



Now we study the convergence of the lefthand sides of ( |2.17|) and ( 2.18 ) as e ^ 0. 
First we show that and Uj, — > |u strong in as e ^ 0. Because of the spherical 

symmetry, compact support and unit integral of ip^ , 

[ dH ^'{l)l®i=\l, (2.21) 

for sufficiently small e. From this result and the definition of it follows that 

uf^(x,t) - -u(x,t) = / d^f (/:)^(£)€(g)^[u(x + £,t) -u(x,t)]. (2.22) 

3 Jj3 

Thus, 

\Wl-\Al'^<I dH^'mvi{- + t)-vi\\L,. (2.23) 

Because u G L^, it follows by a standard approximation argument that ||u(- +t) — ^ as 
^ — > 0. Hence, it follows from (l2.23| ) that ||u|^ — |u||j;^3 — > as e — > 0, as was claimed. Because 
u|, = — u|^, we also deduce that ||u|n — |u||j;^3 — > 0. 
Next we show that — > and pj. — > |p strong in 

L3/2 

as e ^ 0. First, observe that 

p = {—A)^^didj{uiUj), so that p € L^^'^ by the C alder on- Zygmund inequality. From the 
definitions of (pT and pj. it follows easily that 

pf.(x, t) = - • 47r / i'^d£'^%i') ■ -— / d^£ p(x + £, t), (2.24) 

3 Jo \Be'\ Jb^, 

where is the ball of radius £ centered at the origin and [Bil is its volume. Because 93^ has 
unit integral, 

2 2 r°° 1 r 

p^(x, t) - -p(x, t) = - . 47r / /2^£V'(^') • / 'i''^ b(x + ^, t) - p(x, t)], (2.25) 

and thus 

2 2 r°° 1 r 

Wp't - ^p|Il3/2 < - • 47r / i''d£'ip%i') ■—- dH ||p(. +£) -p||^3/2. (2.26) 

3 -3 Jo \Be\ Jb^, 

Consequently, \\pj. — \p\\iz/2 ^ as e — > 0, as was to be proved. Since p\ = p^ — Pj', also 

Wpl - 3pIIl3/2 ^ 0. 



Entirely similar arguments show that ((ux*ux)u)^ — (ux*ux)^u strong in as e — > 0, 
for X = L,T. We leave that to the reader. 

Finally we see that the the lefthand side of equation ( |2.17| ) converges as e ^ to the 
quantity 



(i|u|2+p)u 



in the sense of distributions, and that the lefthand side of ( 2.18 ) converges to | times the same 
quantity in the curly brackets, also in the distribution sense. However, by the result of Duchon- 
Robert ||T[, that quantity in the brackets is equal to minus the defect distribution D{u). We 
conclude then that 

D'k{u)^D{u) (2.27) 

as e — > for both X = L,T, in the sense of distributions. □ 
We now prove the following: 

Corollary 1 Assume that the functions S'l(u, £), S'x'(u, G L"'^([0,T] x T^) defined by 

SL{u,e) = ]- [ diu{i)[6uL{i)f (2.28) 



and 



have limits 



Then, 
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St{u,£) = 1 f du{i)6uL{£)\SuT{£)\^ (2.29) 



Sx(u) =P-lim5x(u,^), X = L,T. (2.30) 



Sl{u) = -^D{u), St{u) = -^D{u). (2.31) 



The proof is quite straightforward. An easy computation for a spherically symmetric test 
function gives 

4 

-Dl{u) = A7T d£ [e^ip'{e)SL{n,ei) + 2£^ip{£)ST{u,ee)] (2.32) 
"J Jo 
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and 

Q rco 

-Df.(u) = 47r / M [^V(^) - 2^V(^)] 5T(u,e^). (2.33) 
Taking the limit as e ^ 0, using Theorem 1, the normahzation of the test function, and the 
above hypotheses, we see that 

^I)(u) = -35i(u) + 25T(u) (2.34) 

and 

^D{vl) = -55t(u). (2.35) 

Solving this linear system gives the result. □ 

In the definition of ^^(u, £) in ( p.29| ) we could instead have replaced l^uj'p by the square 
magnitude of a transverse component 5ut = t*(5u, where t is any unit vector perpendicular to (.. 
Call such a quantity ^^(u,^). However, |(5utP = |(5nTp + [(^n^P, where Su'j- = (tx£)-du. Since 
both transverse components give an equal contribution in the spherical average, S't'(u, £) = 
2S't(u, £). Under the assumptions of the corollary, S'j'(u, ^) then also has a limit distribution 
S't(u) as ^ — > 0, and 

5t(u) = -^^(u). (2.36) 
This is the conventional statement of the "4/15-law". 



3 Discussion 

If we couple our Corollary 1 with Duchon-Robert's characterization of D{u) by the inviscid 
limit, [P, Proposition 4, then we arrive at essentially the following statement: For any "nice" 
spacetime region i? C x [0,T], let 

= lim / / d^xdte(x,t) (3.1) 
•'^o \R\ J Jj^ 

with e(x, t) := z^jVu*^!^ (assuming for simplicity that D{'a!^) = 0). If first ^ and next 
^ ^ 0, then 



{[5uL{t)f)ang,R ~ 'T^R^ (3.2) 
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with {■)ang,R denoting an average of {£, x, t) over 5^ x R. To be precise, the relation should 
hold in the sense of averaging in x,t against smooth test functions 'tp{x,t) rather than against 
a characteristic function l^(x, f). The latter sense is stronger, because, if the result holds for 
characteristic functions of Borel sets R, then, by approximation arguments, it holds also for 
smooth functions. In fact, by the portmanteau theorem Q, it is enough to show that it holds for 
all bounded Borel sets R such that D{dR) = 0, i.e. with no energy dissipation concentrating on 
the boundary. We see from these remarks that the meaning of the Corollary 1 is that the scaling 
observed in the Kolmogorov "4/5-law" H] should hold with just local averaging in spacetime 
and an angular average over the direction of the separation vector. Note that, if the solution 
of the Navier-Stokes equations at finite viscosity is not regular, then the same relation would 
hold with the modified definition e{x,t) := i/|Vu^p + D(u'^), including the dissipation from 
the singularities. 

Another recent work by Nie and Tanveer Q has given a rigorous derivation of both the 4/5- 
and 4/3-laws under conditions similar to those of Duchon-Robert and of the present paper. 
The proof of those laws by Nie and Tanveer is for individual solutions of the Navier-Stokes 
equations (assumed to be strong) without any statistical averaging, and with no assumptions 
of homogeneity or isotropy. In Q , as here, the averaging is over spacetime and over directions 
of the separation vector. However, Nie and Tanveer considered the space average over the 
entire domain and time average over an interval [0, T] with T oo. In contrast, the present 
result holds for any spacetime domain, of arbitrary extent, so long as the Reynolds number is 
sufficiently high, and also for solutions of the Navier-Stokes equations that may be singular. 
In this respect, our result is stronger than that of On the other hand, our proof and that 
in IQ] give no indication how large the Reynolds number must be taken to approach the limit, 
whereas Nie and Tanveer establish precise error bounds for their result. 

There is a slight resemblance of our and Duchon-Robert 's local results with the "refined 
similarity hypothesis" that Kolmogorov postulated in 1962 Q. Kolmogorov considered averages 
eb of dissipation over balls B CT"^ and the structure functions obtained by averaging velocity 
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increments 6u(x,t;£) conditioned upon the local dissipation averaged over a ball -B^(x) of 
diameter £ centered at the midpoint i(x + £). He postulated that these might be calculated by 
his original 1941 theory, e.g. for longitudinal differences. Hence, he obtained for unconditioned 
ensemble averages 

{[5uUi)r) - Cr.{el,)r (3.3) 



Our scaling relation (|3.2| ) is supposed to hold pathwise, for individual realizations. Hence, if it 
holds — along with a uniform and integrable bound for i — then one can infer that 

The similarity is apparent. Aside from the fact that we must integrate over time as well as 



space, the major and significant difference is that in the relation (3.4) the region R must be 
fixed as ^ ^ 0, while in Kolmogorov's hypothesis the ball Bi itself shrinks to zero in the limit. 
Our results and that of Duchon-Robert say nothing about intermittency. 

Although we have proved a rigorous theorem, it is only established under various hypothe- 
ses. We regard those as plausible, but it is nevertheless of some interest to inquire about the 
feasibility of numerical and/or experimental tests of the local results proved here and in In 
1^], the 4/5-law has been verified both in a numerical simulation via a space-average over the 
domain and also in experimental hot-wire data. However, in the same simulation, the 4/15-law 
is not satisfied (S. Chen, private communication). Similar violations of the 4/15-law were ob- 
served in experimental X-wire data from an atmospheric boundary layer [|lO|. In both cases, 
anisotropy of the flow may be suspected of vitiating the result. In that case, angle-averaging 
as employed in the rigorous theorems may improve the agreement. On the other hand, a more 
recent numerical study [|ll| with a 1024^ simulation has found that the regimes of validity of 
both the 4/5- and 4/3-laws, in even their global form, are established only very slowly as the 
Reynolds number is increased. Such a slow approach to asymptotia makes a direct test of the 
local results, especially a verification of the numerical prefactor, rather more difficult. 
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